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Thermalization and collective flow of charm (c) and bottom (6) quarks in ultra-relativistic heavy- 
ion collisions are evaluated based on elastic parton rescattering in an expanding quark-gluon plasma 
(QGP). We show that resonant interactions in a strongly interacting QGP (sQGP), as well as parton 
coalescence, can play an essential role in the interpretation of recent data from the Relativistic Heavy- 
Ion Collider (RHIC), and thus illuminate the nature of the sQGP and its hadronization. Our main 
assumption, motivated by recent findings in lattice Quantum Chromodynamics, is the existence of 
D- and B-meson states in the sQGP, providing resonant cross sections for heavy quarks. Pertinent 
drag and diffusion coefficients are implemented into a relativistic Langevin simulation to compute 
transverse-momentum spectra and azimuthal asymmetries (112) of b- and c-quarks in Au-Au collisions 
at RHIC. After hadronization into D- and B-mesons using quark coalescence and fragmentation, 
associated electron-decay spectra and V2 are compared to recent RHIC data. Our results suggest a 
reevaluation of radiative and elastic quark energy-loss mechanisms in the sQGP. 

PACS numbers: 12.38.Mh,24.85.+p,25.75.Nq 
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Introduction. Recent experimental findings at the Rel- 
ativistic Heavy-Ion Collider (RHIC) have given intrigu- 
ing evidence for the production of matter at unprece- 
dented (energy-) densities with surprisingly large collec- 
tivity and opacity, as reflected by (approximately) hydro- 
dynamic behavior at low transverse momentum (pt) and 
a substantial suppression of particles with high px- This 
has led to the notion of a "strongly interacting Quark- 
Gluon Plasma" (sQGP), whose microscopic properties, 
however, remain under intense debate thus far. 

Heavy quarks (HQs) are particularly valuable probes 
of the medium created in heavy-ion reactions, as one ex- 
pects their production to be restricted to the primordial 
stages. Recent calculations of radiative gluon energy-loss 
of charm (c) quarks traversing a QGP in central Au-Au 
collisions at RHIC have found nuclear suppression factors 
Raa— 0-^-0 A 0,0, comparable to the observed suppres- 
sion of light hadrons at high px, and in line with pre- 
liminary (non-photonic) single-electron (e^) decay spec- 
tra 0. [j. The latter also exhibit a surprisingly large 
azimuthal asymmetry (1)2) 0-001 m semicentral Au-Au 
which cannot be reconciled with radiative energy loss, 
especially if c-quarks are hadronized into Z?-mcsons via 
fragmentation. While the underlying transport coeffi- 
cients exceed their predicted values from perturbative 
Quantum Chromodynamics (pQCD) by at least a factor 
of ^5 |3, energy loss due to elastic scattering parametri- 
cally dominates toward low pt (by a factor \j ^fal H3)- 
But elastic pQCD cross sections [ill ll^ | also have to 
be upscaled substantially to obtain c-quark V2 and Raa 
reminiscent to preliminary data, as shown in a recent 
Langevin simulation for RHIC llOl . In addition, con- 
tributions of bottom (b) quarks |l8l llij will reduce the 
effects in the electron- Raa and -V2 at high pt- Quark 
coalescence approaches suggest that an e ± -W2 in excess 



of 10% can only be obtained if yjj( (i) light quarks im- 
part their V2 on D-mesons (see also Refs. ^a^3), (ii) 
the c-quark V2 is comparable to that of light quarks. 

We are thus confronted with marked discrepancies be- 
tween pQCD energy-loss calculations and semileptonic 
heavy-quark (HQ) observables at RHIC. The resolution 
of this issue is central to the understanding of HQ interac- 
tions in the QGP in particular, and to the interpretation 
of energy loss in general. HQ rescattering also has direct 
impact on other key observables such as heavy quarko- 
nium production (facilitating regeneration) and dilcpton 
spectra (where cc decays compete with thermal QGP ra- 
diation). 

In this letter we investigate the HQ energy-loss prob- 
lem by introducing resonant HQ interactions into a 
Langevin simulation of an expanding QGP. Our calcula- 
tions implement a combined coalescence+fragmcntation 
approach for hadronization, as well as bottom contribu- 
tions, to allow for a quantitative evaluation of pertinent 
e -spectra (t>2 and Raa) which is mandatory for a proper 
interpretation of recent RHIC data. Our main assump- 
tion of resonant D- and B-like states in the sQGP has 
been shown [l^ to reduce HQ thermalization times by a 
factor of ~3 compared to pQCD scattering. Theoretical 
evidence for resonances in the sQGP derives from com- 
putations of heavy and l igh t meson correlators within 
lattice QCD (IQCD) [HEj, as well as applications of 
IQCD-based heavy-quark potentials within effective mo- 
dels HI IH |H Q. Except for the mass and width of 
these states, no further free parameters enter our descrip- 
tion, with degeneracies based on chiral and HQ symme- 
try. 

Heavy- Quark Interactions in the QGP. Following 
Ref. pjj our description of HQ interactions in the QGP 
focuses on elastic scattering, mediated by resonance ex- 
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citations on light antiquarks (q) as well as (nonrcsonant) 
leading order pQCD processes dominated by t-channcl 
gluon exchange. The latter correspond to Born dia- 
grams |25l regularized by a gluon-screening mass m g =gT 
with a strong coupling constant, a s =g 2 /(47r)=0.4. The 
key assumption [18j is that a QGP at moderate temper- 
atures T<2T C sustains strong correlations in the lowest- 
lying color-neutral D- and B-meson channels. Sup- 
port for the relevance of such interactions stems from 
quenched 1QCD computations of euclidean mesonic cor- 
relation functions, which, after transformation into the 
timclike regime, exhibit resonance structures for both 
(heavy) Q-Q and (light) q-q states llfj . In addition, 
applications of 1Q CD-based Q-Q potentials have revealed 
both bound 0, E3 and resonance states jl^l with disso- 
lution temperatures of ~ 2T C , quite compatible with the 
disappearance of the peak structures in the 1QCD spec- 
tral functions. Here, we do not attempt a microscopic 
description of these correlations but cast them into an 
effective lagrangian with vertices ($=D, B), at 

the price of 2 free parameters: the masses of the meson- 
fields, fixed at m£)(s)=2(5) GeV, i.e., 0.5 GeV above the 
Q-q threshold (with quark masses m c ({,)=l. 5(4.5) GeV, 
m u _d—0), and their width, T, obtained from the one-loop 
$ self-energy (which, in turn, dresses the $-propagator) 
with the pertinent coupling constant varied to cover 
a range suggested by effective quark models [2(| 
r=0.4-0.75 GeV. The multiplicity of $ states follows from 
chiral and HQ symmetries alone, implying degenerate 
J p =0 ± and states. We emphasize again that, be- 
sides the mass and width of the $ states, no other free 
parameters (or scale factors) are introduced. 

The matrix elements for resonant and pQCD scattering 
arc employed to calculate drag and diffusion coefficients 
of HQs in a Fokker-Planck approach • Resonances re- 
duce the thermalization times for both c- and &-quarks by 
a factor of ~3 compared to pQCD scattering alone [J^ ■ 

Langevin Simulation. To evaluate thermalization and 
collective flow of HQs in Au-Au collisions we perform 
relativistic Langevin simulations [loj embedded into an 
expanding QGP fireball. In the local rest frame of the 
bulk matter, the change in position (x) and momentum 
(p) of c- and &-quarks during a time step St is defined by 

5x = ^5t, Sp= -A{t,p+5p) pdt+5W(t,p+Sp) (1) 
E 

(E: HQ energy), where SW represents a random force 
which is distributed according to Gaussian noise |34j . 
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The drag coefficient (inverse relaxation time), A, and the 
inverse of the diffusion-coefficient matrix, 
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arc given by the microscopic model of Ref. i includ- 
ing p- and T-dcpendencies (the latter converts into a 
time dependence using the fireball model described be- 
low). The longitudinal diffusion coefficient is set in ac- 
cordance with Einstein's dissipation-fluctuation relation 
to 01 B\=TEA, to ensure the proper thermal equilib- 
rium limit. The latter also requires care in the realiza- 
tion of the stochastic process in Eq. 0); we here use 
the so-called Hanggi-Klimontovich realization j^ij which 
approaches a relativistic Maxwell distribution if the Ein- 
stein relation is satisfied (in the Ito realization, e.g., an 
extra term has to be introduced in A |35|). Finally, 
the HQ momenta are Lorentz-boosted to the laboratory 
frame with the velocity of the bulk matter at the actual 
position of the quark, as determined by the fireball flow 
profile (see below). 

The time evolution of Au-Au collisions is modeled by 
an isentropically expanding, isotropic QGP fireball with 
a total entropy fixed to reproduce measured particle mul- 
tiplicities at hadro-chemical freezeout which we assume 
to coincide with the phase transition at T c =180 MeV |2*gj . 
The temperature at each instant of time is extracted 
from an ideal QGP equation of state with an effective 
flavor degeneracy of iV/=2.5. Radial and elliptic flow of 
the bulk matter are parameterized to closely resemble 
the time-dependence found in hydrodynamical calcula- 
tions [2^ , assuming a flow profile rising linearly with the 
radius. We focus on semiccntral collisions at impact pa- 
rameter 6=7 fm, with an initial spatial eccentricity of 0.6 
and a formation time of to=0.33 fm/c, translating into 
an initial temperature of To =340 MeV. The evolution 
is terminated at the end of the QGP-hadron gas mixed 
phase (constructed via standard entropy balance (28|) af- 
ter about 5 fm/c, at which point the surface flow veloc- 
ity and momentum anisotropy have reached w^=0.5c and 
i'2=5.5% (variations in tq by a factor of two affect the c- 
quark V2 and Raa by 10-20% (less for the spectra), 
while a reduction in the critical temperature to 170 MeV 
increases (decreases) v% (Raa) somewhat less. D-meson 
rescattering in the hadronic phase is neglected). 

To specify initial HQ pr-distributions, Pi n i(pr), and 
especially the relative magnitude of c- and 6-quark spec- 
tra (essential for the evaluation of spectra), we pro- 
ceed as follows: we first use modified PYTHIA c-quark 
spectra with (5-function fragmentation to fit D and D* 
spectra in d- Au collisions [31| . These spectra are decayed 
to single-e ± which saturate pertinent data from p-p and 
d-Au up to p|i~3.5 GeV [Hi EH- The missing yield at 
higher px is attributed to contributions from i?-mesons, 
resulting in a cross section ratio of a b i/ a cS —5-10~ 3 (which 
is slightly smaller than for pQCD predictions and im- 
plies a crossing of c- and 6-decay electrons at pt— 5 GeV, 
as compared to 4 GeV in pQCD). 

Fig. ^ summarizes the output of the Langevin 
simulations for the HQ nuclear modification factor, 
R-AA=Pfin{PT) / Piniipr) {P&n- final ^-distributions) , 
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FIG. 1: Nuclear modification factor (upper panel) and elliptic 
flow (lower panel) of charm and bottom quarks in 6=7 fm Au- 
Au( v / s=200 GeV) collisions based on elastic rescattering in 
the QGP. Red (green) and blue lines are for c- (6-) quarks 
with and without resonance rescattering, respectively, where 
the bands encompass resonance widths of r=0.4-0.75 GeV. 



FIG. 2: Raa (upper panel) and V2 (lower panel) of semilep- 
tonic D- and B-meson decay-electrons in 6=7 fm Au- 
Au(v / s=200 GeV) using HQ coalescence supplemented by 
fragmentation (line identification as in FigEJ. The data ~ 
0,0 are for minimum bias collisions, the Raa from STAR 
for 10-40% central collisions. 



and elliptic flow, V2={{pi — p y )/{pt + Py) ) (evaluated 
at fixed pr)- For c-quarks and pQCD scattering only, 
our results are in fair agreement with those of Ref. |lOj 
(recall that the Debye mass in our calculations is given 
by m g =gT while in Ref. |lOl| it was set to m g =1.5T in- 
dependent of a s ). However, both Raa and vi exhibit 
substantial sensitivity to the inclusion of resonance con- 
tributions, increasing the effects of pQCD scattering by a 
factor of ~3-5. Also note the development of the plateau 
in i> 2 (pt>3 GeV) characteristic for incomplete thermal- 
ization of HQs in the bulk matter. 

Hadronization and Single-Electron Spectra. Semi- 
lcptonic singlc-e* spectra arc a valuable tool to inves- 
tigate heavy-meson spectra in ultrarclativistic heavy-ion 
collisions, since their decay kinematics largely conserves 
the spectral properties of the parent particles 0, 13r3 |. 
To compare our results to measured single-e* in Au-Au 
collisions, the above HQ spectra have to be hadronized. 
To this end we employ the coalescence approach of 
Ref. [ill based on earlier constructed light-quark spec- 
tra |37l |. Quark coalescence has recently enjoyed con- 
siderable success in describing, e.g., the "partonic scal- 
ing" of elli ptic flow and the large p/ir ratio in Au-Au at 
RHIC [37I l38l l3^| . as well as flavor asymmetries in D- 
meson production in elementary hadronic collisions |4fj|| . 



Whereas at low pt most of the HQs coalesce into D- 
and i?-mcsons, this is no longer the case at higher px 
where the phase space density of light quarks rapidly 
decreases. Therefore, to conserve HQ number in the 
B- and Z?-mcson spectra, the remaining c- and &-quarks 
are hadronized using (5-function fragmentation. Finally, 
singlc-e^ Pt- and ^-spectra are computed via B- and 
D-meson 3-body decays, and compared to experiment in 
Fig. |21 We find that the effects of resonances are essen- 
tial in improving the agreement with data, both in terms 
of lowering the Raa and increasing V2- The B- meson 
contribution reflects itself by limiting Raa and V2 to val- 
ues above 0.4 and below 10%, respectively, as well as the 
reduction of V2 above pr— 3 GeV. 

To better illustrate the effects of coalescence we plot 
in Fig. |21 calculations where all HQs are fragmented into 
D- and B-mesons. While Raa is significantly reduced, 
most notably in the Pt— 1-2 GeV region, v 2 also decreases 
reaching at most 6%, which is not favored by current 
data. It is, however, conceivable that modifications in the 
fraction of coalescence to fragmentation contributions, as 
well as improvements in our schematic ((5-function) treat- 
ment of fragmentation, will be necessary once more ac- 
curate experimental data become available. Additional 
corrections may also arise from a more precise determi- 
nation of the b/c ratio and nuclear shadowing. 
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FIG. 3: Same as Fig. H but using HQ fragmentation only. 

Conclusion. We have investigated thcrmalization and 
collective flow of c- and 6-quarks within a relativistic 
Langcvin approach employing elastic scattering in an ex- 
panding QGP fireball in semicentral Au-Au collisions at 
RHIC. Underlying drag and diffusion coefficients were 
evaluated assuming resonant D- and B-meson correla- 
tions in the sQGP, enhancing heavy-quark rcscattcring. 
Corresponding pT-spcctra and elliptic flow of c-quarks 
exhibit a large sensitivity to the resonance effects, low- 
ering Raa down to 0.2 and raising V2 up to 10%, while 
the impact on 6-quarks is small. Heavy-light quark co- 
alescence in subsequent hadronization significantly am- 
plifies the V2 in single-electron decay spectra, but also 
increases their Raa, especially in the pr— 2 GeV region. 
Bottom contributions dominate above 3.5 GeV reducing 
both suppression and elliptic flow. The combined effects 
of coalescence and resonant heavy-quark interactions are 
essential in generating a uf of up to 10%, together with 
R AA ~0.5, supplying a viable explanation of current elec- 
tron data at RHIC without introducing extra scale fac- 
tors. Our analysis thus suggests that elastic rescattering 
of heavy quarks in the sQGP is an important compo- 
nent for the understanding of heavy-flavor and single- 
electron observables in heavy-ion reactions at collider en- 
ergies. While induced gluon-radiation is expected to be 
the prevalent interaction with the medium at sufficiently 
highpT QUI [Ol' it may not be the dominant effect below 
Pt— 6 GeV or so. A complete picture should clearly in- 
clude both elastic and inelastic rescattering mechanisms. 
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